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Abstract 

We show that a 2N junction SQUID (Superconducting QUantum Interference 
Device) made of 2N overdamped, shunted, identical junctions may be described as 
a system having only 6 degrees of freedom for any N > 3. This is achieved by 
means of the reduction introduced by Watanabe and Strogatz (Physica D, 74, (1994) 
197) for series biased arrays. In our case 6 rather than 3 degrees of freedom are 
necessary to describe the system, due to the requirement of phase quantization along 
the superconducting loop constituting the device. Generalization to multijunction 
parallel arrays is straightforward. 



1 Introduction 

Two linear arrays, each containing N Josephson junctions, closed into a superconduc- 
tor loop constitute a 2N-junction SQUID (Superconducting QUantum Interference De- 
vice), (see Fig.l). Such multijunction structures are of interest in designing high gain 
SQUIDs@, il|. The case with a single junction for each branch (N=l) is the well 
studied case of a dc-SQUID[j5], |6, |7J. Structures containing N=2, and N=3 junctions 
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for each branch, have also been fabricated and investigated in some detail Jl], . More- 
over theoretical analysis based on numerical simulations of the 2N-junction SQUID has 
appeared toop]. The aim of this letter is to point out that, as for series arrays of over- 
damped Josephson junctions shunted by passive elements, a 2N-junction SQUID, made of 
overdamped (vanishing capacitance) identical junctions, exhibits a highly degenerate dy- 
namics. In particular we will show that the transformation introduced by Watanabe and 
Strogatz (WS) |8] may be fruitfully applied to the case considered here, with few modifica- 
tions, to demonstrate the existence of 2N — 6 constants of motion, i.e., that any trajectory 
in the phase space of a 2N SQUID is confined to a 6-dimensional subspace. 



2 The 2N-junction SQUID model 



The system we are interested in is sketched in Fig. 1 as a superconducting loop in which 
each of the two branches include N identical Josephson junctions. The behavior of the 
circuit is described by 2N dynamical variables, {0^ }, {<$'} (& = h ---N) representing the 
gauge invariant phase difference across each junction for the branch (a) and (6), respec- 
tively. In the case of overdamped junctions, the governing equations are 000 
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R is the normal resistance, J 0 is the critical current, / is the total bias current. The coupling 
terms J and — J coincides with the screening current, which is determined through the 
fluxoid quantization: 
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Here $o is the flux quantum, 2 L is the inductance of the loop, and <f> a is the externally 
applied magnetic flux. The key observation is that the phases 0^ and are globally 
coupled through Eq. (3) and that for each subsystem, (a) and ( b ), this coupling is inde- 
pendent of the index k. This generalizes the case of globally coupled junction series arrays 
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studied by WS [;8j. Indeed, although our system is not globally coupled (there is a change 
of sign of the coupling term when passing from a branch to another), it can be split into 
two globally coupled sub-systems . Each subsystem contains N series connected junctions 
and satisfies the conditions for the applicability of the reduction introduced by WS. If we 
introduce the dimensionless variables 
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Eq.s (1-3) take the form (dots denote derivative with respect to the normalized time): 
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Therefore we introduce a change of coordinates separately for each branch, as 
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Here ©( a ) ; y( a ) ) vld a ),0 ( H, and vj/( b ) ; are unknown functions of time (0 < 7^ , 7® < 
1), and are 2N constants. It is a straightforward exercise to prove that the 
transformations (7,8) reduce the 2N Eq.s (4-6), for every N > 3, to the following six 
equations: 



7 (a) = -(l- 7 ( Q)2 )cos0 (a) , 



(9) 
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7 (a) ^ (a) = yjl - 7 W 2 sin 0 (a) , ( 10 ) 

7 (a)@(a) = 7 (o) Q~j) +sm 0 (a) , ( 11 ) 

¥ b) = -(i-7 (fe)2 ) cosO (b \ (12) 

7 (b) T (6) = yj\ -7W2 sin 0(6) > (13) 

7 ( 6 )@( 6 ) — 7^ ^ + sin 0^. (14) 



Eq.s (9-14) form a closed system for the six unknown functions, since j can be expressed 
in terms of the new variables inserting the change of coordinates (7,8) in Eq. (6). The 
’frozen phases’ (adopting the WS’s language) appear in the system merely as parameters. 
Some comments are in order. Firstly, the above system becomes singular if either 7 ^ a \ or 
7® vanish. Moreover negative values for 7^, or 7^ seem not excluded, contrary to their 
definition. WS have shown that these troubles are not essential, but just an artifact of 
the chosen coordinate system. As a matter of fact they can be avoided by the additional 
change of coordinates: 
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Thus the system, in terms of the new variables, becomes : 
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(25) 
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which is well-behavied even at 7 ^ = 0, or 7 ^ = 0 The second point is about the initial 
conditions. Given 2N arbitrary initial junction phase values, the values for the 2N ’frozen 
phases’ and the initial values for the six ’reduced variables’ x^^ b \ z ( a ),( b ) 

should be obtained from Eqs. (7,8). This of course can be done in many ways: a simple 
and natural choice is the ’identity transformation’ : 



X G)(0) = 2 /(°>( 0 ) = z (a) (0) = x^ b \0) = y {b) ( 0) = z (6) (0) = 0 



(27) 



which straightforwardly gives : 

=<fi am (0) (28) 

The investigation of the meaning of such moltiplicity of choice for the initial values is an 
important issue for the understanding of the dynamics of the system, but goes outside the 
scope of this letter. We just comment that, as (27) and (28) show, there is at least one way 
to effectively achieve the reduction. 

We have numerically verified, for few values of the parameters, that the reduced system 
of six equations yelds the same results of the originary system, Eqs. (4-6). This was checked 
with increasing number of junctions up to N= 20 junctions per branch. Care must be taken 
in reconstructing the phase values from the six unknown functions to avoid spurious 2tt 
jumps in the solution. 



3 Parallel arrays 



The above property may be exploited, more generally, to a multijunction parallel array of 
M branches each branch containing N identical junctions (Fig. 2). Then Eq.s (1-2) can 
be generalized as follows: 
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where / ^ is the current in the m — th branch, and is the phase of the k-th junction 
in the m-th branch. The I ^ sum, of course, to the total bias current I : 
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Moreover the condition for the fluxoid quantization in each loop writes : 
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Here c\^ is the mutual inductance coefficient accounting for the contribution to the mag- 
netic flux at the loop (j, j+l)due to the current flowing in the l - th branch. These coefficients 
are determined by the specific geometry of the circuit. Eqs. (30) and (31) are a set of 
linear equations for the unknown I^ m \ and allow a unique determination of their value as 
function of all the phase variables. In the simple case in which only the nearest neighbours 
branch currents contribute to the magnetic flux in each loop m? 0- cJ j) assume the simple 
form : 



(31) 
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where L is the inductance of each branch and is the Kronecker’s symbol. Under these 
simplifying hypotheses it is easy to show that : 
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The N x M Eqs. (29) and (33) describe completely a multi-junction parallel array, 
and, as in the case of the 2Nj-SQUID, it can be viewed as M sub-systems each one globally 
coupled to the rest of the system. The voltage V developed across the circuit may be 
obtained by time derivating Eq.(33): 
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It should be pointed out that the simplification introduced here (to neglect off diagonal 
terms of the mutual inductance matrix) is not essential to obtain a global coupling, i. e. also 
the more general case will lead qualitatively to the same kind of (global) coupling. As in 
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the previous discussion of the 2N-junction SQUID, it is therefore possible to introduce M 
transformations analogous to (7,8) (one for each superconducting branch) and to conclude 
that only the following 3 M equations (m = 1...M) 
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should actually be solved for the 3 M variables { x ^ m \ y ( m \ z rather than N x M. It 
should finally remarked that the all identical junction assumption may be partially relaxed. 
Indeed no particular troubles are met in the previous derivation if one assumes the junctions 
in a branch (or even their number) to be different from those of another branch. That is 
to say, the requirement of identical parameters concerns only the junctions in the same 
branch. 



4 Conclusion 

The possibility to reduce an N x M system to 3 x M one is in itself of some advantage 
to simplify numerical simulations. Still we believe that the major achievement that we 
have reached is that the low dimensional motion for ’globally coupled’ systems with sine 
nonlinearity can be extended also to sets of equations that are not ’globally coupled’ in a 
strict sense, but can be divided in subsystems exhibiting global coupling (at the price to 
increase the number of dynamical variables necessary to describe the system). So we have 
explicitly proved that phase quantization in superconducting loops has no other effect that 
to introduce 3 more variables for each loop. Beside the extension of the theory there is 
also a consequence of practical importance: neutral stability is a drawback for practical 
applications (as local oscillators) because of the consequent deterioration of the lincwidth 
of the emitted microwave. Our findings can be interpreted as a criterion for the maximum 
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number of junctions that can be inserted in a loop: above 3 the system undergoes a 
reduction similar to that illustrated in this letter. It should be noted that this criterion 
is not a positive one, i.e., there might be configuration that are neutrally stable but do 
not contain more than 3 junctions per loop [11, 1^]. At this stage it is not known to 
which extent the formal analogy between the system studied in Ref. [|8j] may be carried 
on. Indeed, there are differences in the actual form of the coupling term that will lead to 
substantial differences: for instance Eq.s (4-6) do not exhibit the important property of 
reversibility (changing in —(j)^ and t! in —t' does not leave unchanged the equations). 
This will lead to profound effects on the dynamics of the system. 
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Figure Captions 

Fig.l Schematic circuit model for a 2N-junction SQLIID; 

Fig. 2 Schematic circuit model for a multi- junction parallel array; 
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